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Abstract 

Consider a continuous time Markov chain with rates Q in the state space Au{0}. 
In the associated Fleming- Viot process N particles evolve independently in A with 
rates Q until one of them attempts to jump to state 0. At this moment the particle 
comes back to A instantaneously, by jumping to one of the positions of the other 
particles, chosen uniformly at random. When A is finite, we show that the empirical 
distribution of the particles at a fixed time converges as N — > oo to the distribution of 
a single particle at the same time conditioned on not touching {0}. Furthermore, the 
empirical profile of the unique invariant measure for the Fleming- Viot process with 
N particles converges as N — > oo to the unique quasi-stationary distribution of the 
one-particle motion. A key element of the approach is to show that the two-particle 
correlations is of order 1/N. 
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1 Introduction 

Let A be a finite or countable state space, whose elements are also called sites. Let 
Q = (q( x jU)] x,y ^ AU {0}) be the transition rates matrix of an irreducible continuous 
time Markov jump process on A U {0}. 
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For each integer iV > 1, the Fleming- Viot (fv) process with iV particles, is a continuous 
time Markov process, £ t e i > 0, where denotes the position of the particle z at 
time t. The generator C N acts on functions / : — > R as follows 

TV 

(i-i) 

i=l X £A\{((i)} 

where i %,x {i) = x, and for j ^ i, £ l ' x (j) = £0)- We set [/(£*)] = exp(tC N )f(^). In words, 
each particle moves independently of the others as a continuous time Markov process with 
rates Q, but when it attemps to jump to state 0, it comes back immediately to A by 
jumping to the position of one of the other particles chosen uniformly at random. 

This type of FV process was introduced by Burdzy, Holyst and March in [1] for Brownian 
motions on a bounded domain. In their model N brownian particles evolve independently 
until one of them reaches the boundary, which plays the role of state 0. In [I], the authors 
prove that the empirical profile of the invariant measure converges in this case to the 
first eigenfunction of the Laplacian on the domain with homogeneous Dirichlet boundary 
conditions. 

Denote by ??(£, x) the number of £ particles at site x, and by m(£) the empirical measure 
induced by a configuration £ G UatA^ 

i^,*) :=£>{£(•)=*} ^d m(0:=%^%^- (1.2) 

We also use m x (£) to denote m(£)(x) and q(x, x) = —J2 y eAu{o}\{x}l( x ^y)- With this 
notation, the time-derivative of E^m^^t)] is easily seen to be 

Consider the process on A U {0} generated by Q, with initial law \i and denote T t [i its 
law at time t conditioned on not having touched {0} up to time t. In other words, for all 

x e A 

(x) = E, € A^(y)exp(tg)(y,x) 

1 - E^eA exp(tQ) (y, 0) ' 
Let .M be the space of probability measures on A. Then {T t ,t > 0} is a semi-group on 
M. and T t fi is the unique solution to the Kolmogorov forward equations: T (/i)(x) = fJ,(x), 
and 

■jT t n{x) = q{v, x ) T t fi(y) + Y q ( y ' °) ^Wy) T t/ i ( x )- ( L5 ) 

ygA j/eA 

A quasi- stationary distribution (qsd) for Q is a probability measure v on A that is 
invariant under {T t ,t > 0}, that is 

Tw = za for all t > 0. 
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Now, assume that [i is close to m(£), and look at (jl.3p and (11.51) . A natural approach 
to show that T t fi is close to Ef [m(£i)] is by establishing that the occupation numbers of 
two distinct sites, at time t, become independent when N tends to infinity (the so-called 
propagation of chaos). For this purpose, Ferrari and Marie [3], estimate the correlation of 
two ^-particles, when A is only assumed countable. 

Proposition 1.1 (Ferrari and Marie [3], Proposition 3.1). Let \x be any probability measure 
on k, and /i® N the product probability on A N . Then, there is a constant k, such that for 
any x, y G A, and t > 0, 



Ef K(6)m y (6)]^(0 - J EfK(6)]rf^(0 J Ef[m y &)]dii® N (0 



e 

< — 

- N 



;i.6) 



For countable A, the ergodicity of the FV process is not guaranteed, neither is the 
convergence of the invariant measures as N — > oo. Under further assumptions (see (11.81) ). 
there exists a unique invariant measure for the FV process and its profile converges as 
N — > oo to the unique QSD. 

Theorem 1.1 (Ferrari and Marie [3], Theorem 1.2 and 1.4). Let fi be any probability 
measure on A, and fi® N the product probability on A N . Then, for each x G A 

lim / Ef [K(6) - T tf i(x)f] d^ N (C) = 0. (1.7) 

Moreover, if 

\ inf q(x, z) > maxg(x,0), (1.8) 

1 xGA\{z} x£A 

then there exists a measure v on A and for each N there exists a unique stationary measure 
for the Fleming- Viot process, say X N such that for each x G A 

lim / (m x (0-v(x)) 2 d\ N (O = 0- (1-9) 
N ~*°° Jan 

Furthermore, v is the unique quasi-stationary distribution of Q. 

Here, we have similar results. We prove the asymptotic independence of two particles 
but, in contrast to [3], we consider deterministic initial configurations and obtain bounds 
on the correlations that hold uniformly on the initial distribution of the particles. This 
result also holds for countable A. 

Proposition 1.2. For each t > 0, and any x, y G A 

sup |EfK(6K(6)]-Ef[m x .(6)]EfK(6)]| < %^ + ^(e 2Ct -l). (1.10) 



Our proof yields also that any finite number of particles evolve independently in the 
limit N — > oo. A similar result is obtained by Grigorescu and Kang [5], following the 
approach [5] by the same authors. 
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Finite state spaces. In the rest of the paper we consider a finite A. In this case, for 
each N > 2 the FV process is an irreducible pure-jump Markov process on the finite state 
space A^. Hence it is ergodic, that is, there exists a unique stationary measure for the 
process and starting from any measure, the process converges to the stationary measure. 
We still denote this measure with A . 

When A is finite, Darroch and Seneta [2] prove that T t fi converges exponentially fast 
to a probability measure u, uniformly in the initial measure. The measure v is the unique 
quasi- stationary distribution of Q. 

Theorem 1.2 (Darroch and Seneta, 1967). Assume A is finite, and that the process on A 
with rates {q(x,y),x,y G A} is irreducible. Then there exists 9 > 0, such that 

sup \\TtH - v\\ < e~ 9t . (1.11) 

We have used the total variation norm \\fj, — u\\ = J2 xeA — u ( x )\- 

The asymptotic independence of Proposition 11.21 implies naturally the convergence of 
the empirical means in the FV process to the conditioned distribution T t m(£), uniformly 
in £. Moreover, A finite yields that Tf(m(£)) is close to the unique QSD, uniformly in £ as 
implied by (11.111) . These two facts imply the following result. 

Theorem 1.3. Assume A finite. For any positive time t 



lim sup Ei v [||m(&) -T t m(^)\\) = 0. (1.12) 



Moreover, 



N^oo 



lim / \\m{£)-v\\d\ N {£) = 0. (1.13) 



A iV 



Remark 1.1. Note that the limit f II . 1 3 j) follows from (11.121) using (11.111) . 
1^(0-^11^(0= / E?[\\m(t t )-v\\]d\ N {t) 



< 



[ Ef[\\m^ t )-T t m(0\\}dX N (0+ [ \\T t m{t) - u\\ d\» (£) (1.14) 
Ja n Ja n 



< sup E € [||m(&) - T t m(0l|] + sup \\T t /i - u\\. 

£eA N fj.eM 

On the other hand, (I1.13P implies readily (using that < m x < 1 for all x G A), that for 
any subset U C A 

lim !X\m x d\ N = T\v{x). (1.15) 

J x&j xeu 
We omit the obvious inductive proof. 

The paper is organized as follows. In Section [2], we construct the process a la Harris fol- 
lowing [5] . We use this construction to estimate the correlations and prove Proposition 11.21 
in Section [31 Finally, in Section HJ we prove Theorem II. 3 1 
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2 Graphical construction 



The graphical construction is used to prove the asymptotic independence property. A 
realization of the process t > 0) is a deterministic function of a realization of a marked 
point process. All initial conditions £ use the same realization of the marked point process. 

Let C := max x€ A q(x, 0) be the (maximum) absorption rate, and 



To each particle i, we associate two independent Poisson marked processes (u>{,cuY), which 
we call respectively the internal and voter point processes, described as follows. 

• The internal process is defined on R x A A with intensity measure qdt d r y(F) with 



that is, 7 is the joint distribution of independent random variables with marginal 
distributions {p(x, ■), x G A}, so that to each state x the random mark F assigns 
a state y = F(x) with probability p(x,y). An internal marked-time is {t,F), which 
means that if at at time t— , particle % is at site x, then at time t, particle i jumps to 
site F(x). This gives the correct rate q(x,y) = qp{x,y) for jumps from x to y. 

• The voter point process is defined on R x ({1, . . . , N} \ {i}) x {0, 1} A with measure 
intensity C dtd@i(j) dj'(Q, where $ is the uniform probability on {1, ... , N} \ {i}, 
and 7' is the joint distribution of independent Bernoulli random variables with pa- 
rameters g ^°- > , x G A. A voter marked-time is (t,j, C), where j corresponds to a 
reincarnation label, and ( takes into account the position-dependent rate: if the z-th 
particle is at position x at time t—, it jumps to the position of particle j at time t 
only if ((x) = 1, yielding the correct rate q(x, 0)/(N — 1). 

We call uj = ((a;/ , i e {!)•••) ^}) an i-i-d- sequence of stationary marked point pro- 
cesses associated with labeled particles. Finally, for any subset of labels a C {1, . . . , N}, 
we denote by u a the processes associated with labels a. For any s < t reals, we denote by 
u; a [s,i[ and u; a [s, £] the projections of the marked-times in the time period [s,t[ and [s,t], 
respectively. 

We construct > 0} in such a way that £ f is a function of the initial configuration 
£0 and the time-marks u[0, t], t > 0. Fix an initial configuration £ G A^, and t > 0. There 
is, almost surely, a finite number of time-marks within [0, t], say i^, and let {bk,k < K} 
be the ordered time-realizations. We build £ f inductively as follows. 

• At time 0, the configuration is £o- 



g:=sup ^2 i( x iy)', P( x iV) : 



q{x,y) 



,y^x; p(x,x):=l- ^ P( X 'V)- 

y£A\{x} 




d-f(F) = Y[p(x,F(x)), VFgA a . 
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Assume £b fe is known. For t G [bk, bk+i[, set £ t = £,b k - We describe now & 



'fc+r 



— If bk+i corresponds to an internal time of particle i and mark F, we move particle 
i to F(x) where x = £b fe (z). This move occurs with rate qx q ^ x,F j- x ^ = q(x, F(x)). 

— If bk+i corresponds to a voter time of particle i and mark (J, Q, we move particle 
i to the position of particle j if ((x) = 1 where x = £b k (i)- This move occurs 
with rate C x ^ x = q fa o)/(iV - 1). 

It is easy to check that {£t,t > 0}, as constructed above, has generator given by (11. ip . 

By translation invariance of the law of u, if we use the marks ou[— t, 0] instead of the 
marks t], the configuration so obtained has the same law of £t as constructed above. 
We abuse notation and call £t the configuration constructed with the marks uj[— t, 0]. For 
each particle label i, we build simultaneously a set of labels i^ l {t) of particles which could 
potentially influence £t(i); this set is also a function of u>[—t, 0]. First, the process 1 1— > ip % {t) 
may only change at the time-realizations of the voter process to , and it changes as follows. 
Let — v be the largest time realization of ujj[—t, 0[, and let (j, Q its associated mark. Then, 
for < s < v, we set if) l (s) = {i}, and ij) l {v) = {i,j} (regardless the values of £). For s < t, 
assume that ip % (s) is built, and let — v be the largest time realization of u^ t ^[—t, —s[, and 
let (j, Q its associated mark. Then, 

Vu G [s, u[ = and Vfy) = ^(s) u 0'}- 

Note that for any t > 0, ^ l (t) is o^cj^f— £, 0[)-measurable, and that for any subset of labels 
a C {1, . . . , iV} containing z, we have 

= a} G ff(^[-t, 0[), and {^(t) = a, &(z) = x} G a(a; [-t, 0[). (2.16) 

The next lemma says that the sets of labels associated to two different particles intersect 
with probability of order 1/N. 

Lemma 2.1. For i,j distinct labels, and t > 

P n &(t) + 0) < (e 2Ct - 1) . (2.17) 

Proof. First, we show that the rate of growth of ip l (t) is at most exponential. It is clear 
from the construction of ip 1 , that its rate of growth at time t is at most Cip l (t), and that 
it grows by adding one label (from {1, . . . , N} \ ip l (t)) uniformly at random. Thus, 

1 \ J < c\nt)\ =► nmm < exp (ct) . (2.1s) 

Second, we show that for i,j two distinct labels, 

P (nt) n ^'(t) ^ 0) < ^-y j( E[^( S )]E[^'( S )]d S . (2.19) 
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Note that 

P(^(t) n^'(t) ^ 0) = / E 

and 
E 



d¥ (^(s) n + 0|ff(w[-s, 0[)) 



r/s 



(2.20) 



dP (^(s) n ^ 0|<r(w[-s, o[)) 



2C 



N 



-E [V( s )n^( s )=0l^X s )ll^(> 



2C 



AT 



- \a\\b\¥(^(s)=a,^(s) = b) 



anb=<D 



2C 



N 



- |a||6|P(V'(«)=o)P(^(«) = 6) 



anfe=0 



< 



2C 



'TV 



T E [^(s)] E [^(s)} 



(2.21) 



We used in (I2.2ip that for two non-overlapping subsets of labels a and b, {ip l (t) = a} and 
{ip^(t) = b} are independent by (12.161) . 

This concludes the proof of (12351) . Now, (l2~T7j) follows from (12~T9|) and (1215)1 . ■ 



3 Proof of Proposition 11.2 



We need to show that for any x,y E A, any time t > 0, and initial configuration £, 

|E[r/(6)(x)r / (6)(y)]-E[r / (6)(a;)]E[r / (6)(2/)]| < 2Ne 2Ct . (3.22) 

Here and throughout this section we use E and P to denote E^ and P^ respectively. Using 
(ll.2p . the difference of expectations in the left hand side of (I3.22p is 



i<N j<N 



(3.23) 



For a subset a, = x,ip l (t) = a} is cr(u; a [— t, 0[)-measurable, by remark (12.161) . Thus, 

for two non-overlapping subsets of labels a and b 

P(^(t) = a^{t) = b, = = y) 

= P (^(t) = a,&(i) = s) P (^'(t) = 6,6(j) = y) • 
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Compute a generic term on the right hand side of (I3.23P with i ^ j: 

mm = x,ttU) = y}) = nnt) n + %m = ^mj) = v) 

+ ]T p (^(0 = oi^{t) = b,m = *,ttU) = y) 

aflb=0 

= p(^(t) n v(t) ^ 0,6(0 = x,Uj) = v) 

anb=0 

(3.24) 

To compute P(6(0 = x)P(6(j) = y), we can think of two independent marked-point 
processes driving the evolution (put a tilda for the independent copy), and we have a 
decomposition similar to (13.241) for i ^ j: 

nm = x)nm = y)= m\t) n ± 0,6(0 = xMj) = v) 

+ J2 p (/W = a ^(f) = b M*) = x Mj) = y) 

anb=<6 

= P(^(t) n ft(t) ^ 0,6(0 = x,6(j) = 2/) 

an6=0 

(3.25) 

Subtracting (13.241) and (I3.25P we get that for i ^ j, 

\mm =x,6(j) = y}) - mm = x)nm = y»i 

< P(#)n^(i)^) + P(#)n^(^0) < JL_(e 2Ct -i), 

by Lemma 12.11 (we have used that the Lemma holds also for i[} l (t) D i/ji(t)). Thus, by 
summing over % and j G {1, . . . , N}, and noting that there are N diagonal terms which 
bring a factor N when i = j, we get the desired bound. 



4 Proof of Theorem II. 3L 



In view of (11.141) . we first estimate [||^(6) — ^~t m (0)ll]- ^ i s m ore convenient to work 
with / 2 -norm, rather than total variation norm. For a function tp : A — > R, we denote its 
Z 2 -norm as 

IM| 2 =(£^(*))'- 

By Cauchy-Schwarz, note that if /x, z/ are probabilities on A, 



| yLi — ^|| 2 < Wf 1 ~ HI < V |A| — Z/||2. 
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To estimate Ef [||m(£t) — T t m(£) H2] note that 

Ef [\\m(Ct)-T t m(Oh] < Ef [||m(&) - Ef [m(6)] || 2 ] + ||Ef [m(^)-T t m(C)h- (4-26) 
Taking w = x in f ll.lOp we obtain 

2e 2Ct 

E [K(6) - EKfe)]) 2 ] < - ir . (4.27) 
By (14.271) and Jensen's inequality, we have 

[Ef |K&) - Ef [m(£ f )] || 2 ] 2 < Ef [||m(6) - Ef [m(£ t )] H*] < . (4.28) 

The second term in (14.261) is dealt with in the following lemma. 
Lemma 4.1. For any T > 0, 

Urn max max ||Ef [m(£ t )] - T t m(0h = • (4.29) 

Proof. We introduce some simplifying notations. 

u x (t)=Ef[m x {£ t )], and « x (*) = T t m(£){x). (4.30) 
We show that there is a constant S such that for any t > 

7 4e 2CI * a(v 0) 

-\\u(t) - v(t)\\l < B\\u(t) - v(t)\\l + \ (4.31) 

Since \\u(0) — i; (0) || 2 = 0, the result follows at once by means of Gronwall's inequality. 
We fix t > 0, and we often omit to display the time-dependence. From f jl.3D 

= ^2 q ( y > x "> u y + q ( y ' °) u * u y + R *(£> *) ( 4 - 32 ) 

where t) = ^ ? (y, 0) ^yEf [m, (6)7^(6)] - Ef K(6)] Ef [m x (&)] 

Proposition ll.2l implies that 

sup sup li^^t) I < ■ (4.33) 



xeA £ 

On the other hand, from (I1.5P 



AT 



^2q(y,x)v y + ^2q(y,0)v x Vy. (4.34) 

2/eA yeA 



Subtracting (14.341) from (I4.32p . we obtain 

d(u x - v x ) 



Y Q(v> x ) K - v v) + Y ?(2/» °) ~ Wy) + £) • 

j/GA «/GA 



Now, 



lll wt) _„ (i) l li=E *^) (IiI _„J 

X 

= YY q ( x,y " > ( u y ~ v y> ( Ux ~ Vx "> + YY q ( y ' °) _ W:r ^) ( M:r ~ ^ 

rrGA y<=k x-gA 2/GA 



+ ^i? :c (u :c - v x ). 



(4.35) 



We deal with each term on the right hand side of (14.351) in turns. First, 

1/2 



\YY q ^ x ' y ^ u y~ v y^ Ux ~ Vx ^ - I Y q2 ( x >y) I 11^ — ^lll- ( 4 - 36 ) 

xGA yeA \x,yeA J 

To deal with the second term observe that u x u y — v x v y = v x {u y — v y ) + u y {u x — v x ) and 
hence the second term in (14.351) term equals 

YY q ( y ' Q ) v *( u y ~ v y)i u x ~ v x ) + } j } j q{y, 0)u y (u x - v x )(u x - v x ). (4.37) 
xga yGA xga j/ga 

Since v is bounded, the first term in (14.371) can be treated as (14.361) with q(x, y) replaced 
by q(x, y) := q(y, 0)v x . The second term in (14.371) equals 



Y q ( y ' °"> u y - W U ~ V W 2 2 SU P q ( y ' °)- 



Finally, taking into account that J2 x v x — J2 x u x — 1 ; the last term of (14.351) can be 
bounded by 

o , d , ^ Ae2Ct «to> °) 

2 sup \R X \ < f 

x J' 

Collecting all these computations, we obtain ( 14.311) . ■ 
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